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Quantitative descriptions of network strncture in big data can provide fundamental insights into 
the function of interconnected complex systems. Small-world structure, commonly diagnosed by high 
local clustering yet short average path length between any two nodes, directly enables information 
flow in coupled systems, a key function that can differ across conditions or between groups. However, 
current techniques to quantify small-world structure are dependent on nuisance variables such as 
density and agnostic to critical variables such as the strengths of connections between nodes, thereby 
hampering accurate and comparable assessments of small-world structure in different networks. 
Here, we address both limitations with a novel metric called the Small-World Propensity (SWP). 
In its binary instantiation, the SWP provides an unbiased assessment of small-world structure 
in networks of varying densities. We extend this concept to the case of weighted networks by 
developing (i) a standardized procedure for generating weighted small-world networks, (ii) a weighted 
extension of the SWP, and (iii) a stringent and generalizable method for mapping real-world data 
onto the theoretical model. In applying these techniques to real world brain networks, we uncover 
the surprising fact that the canonical example of a biological small-world network, the C. elegans 
neuronal network, has strikingly low SWP in comparison to other examined brain networks. These 
metrics, models, and maps form a coherent toolbox for the assessment of architectural properties in 
real-world networks and their statistical comparison across conditions. 


I. INTRODUCTION 

In the era of big data sciences, the ability to infer fun¬ 
damental principles of system function from large sets 
of interconnected variables is of increasingly pressing im¬ 
portance for structural characterization, functional pre¬ 
diction, and novel network design. Traditionally, net¬ 
work science has been posited as a particularly appealing 
framework in which to form such inferences based on its 
simplicity and ubiquitous utility [l|. Yet, it is exactly this 
simplicity - arguably network science’s greatest strength 
- which often ignores critical system details, causing it 
to dual as the field’s greatest potential weakness 0 . For 
example, many network-based tools were originally de¬ 
veloped to examine isolated networks constructed from 
binary links. However, storage facilities around the world 
now house rich new data offering multiple instances of the 
same or related networks, with finely measured weights 
on links between nodes. Progress in using these data 
to understand and controllably manipulate complex sys¬ 
tems requires complementary theoretical advances in the 
realism of network-based analysis tools. 

Particularly relevant tools are those used for charac¬ 
terization: quantitative statistics to describe the organi¬ 
zation or structure of networks, from which one can infer 
properties of their function. A quintessential example is 
the structure of small-worldness, ubiquitous across many 
real-world networks Q, whose local clustering combined 
with the ability to move quickly through the network has 
been shown to have important implications for functions 
from synchronizability to information flow Watts 


and Strogatz formalized the concept of a small-world net¬ 
work by describing a simple theoretical model on binary 
networks Q. Nodes are placed on a lattice and connected 
to their nearest neighbors within a radius, r, giving the 
network a high clustering coefficient. Shortcuts are then 
introduced to the network by randomly rewiring each 
edge with some probability, p. For p = 0 the network 
is a lattice, for p = 1 the network is random, and for a 
range of intermediate p, the network is small-world. 

While this model has proven useful for theoretical in¬ 
vestigations dSl , it does not provide a direct metric to 
assess small-worldness in observed real-world networks. 
If one has observations of a network over a range of 
size scales, one can determine small-worldness by ask¬ 
ing whether the characteristic path length scales as the 
loo(iV) where N is the number of nodes in the network 
[1^ . More commonly, however, one has observations of 
networks at a single size scale, and therefore needs to map 
real-world data to the theoretical Watts-Strogatz model. 
Such a mapping can be achieved by a comparison of the 
observed clustering coefficient andpath length to that of 
random and/or lattice networks |14l - ll^ . 

Critically, these models - and associated tools and 
mappings - have lagged behind advances in the data 
sciences enabling more comprehensive network measure¬ 
ment. Specifically, these statistics are dependent on 
network density and neglect critical variables such as 
the strengths of connections between nodes, limiting 
their ability to diagnose and compare small-world struc¬ 
ture in different networks. To address these limita¬ 
tions, we introduce a novel diagnostic called the Small- 
World Propensity (SWP), which quantifies the extent 
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to which a network displays small-world characteristics 
while accounting for variation in network density. We 
next present (i) a standardized procedure for generat¬ 
ing weighted small-world networks, (ii) a weighted exten¬ 
sion of the SWP, and (iii) a stringent and generalizable 
method for mapping real-world data onto the theoretical 
model. We verify our methods using standard benchmark 
networks and use them to examine small-world proper¬ 
ties of real-world brain networks in humans and neuronal 
networks in worms. Surprisingly, we observe that the 
neuronal network of C. elegans, originally touted as a 
quintessential example of a biological small-world net¬ 
work, shows especially weak SWP. 

II. RESULTS 
A. Small-World Propensity 

To quantify the extent to which a network dis¬ 
plays small-world structure, we define the Small-World 
Propensity, (/), to reflect the deviation of a network’s clus¬ 
tering coefficient. Cobs, and characteristic path length. 
Lobs, from both lattice {Ciatt, Liatt) and random {Crand, 
Lrand) netwoi'ks constructed with the same number of 
nodes and the same degree distribution: 


As illustrated by the application to the Watts-Strogatz 
model, the SWP is best utilized as a comparative metric 
to describe the extent to which a network displays small- 
world structure. While a definitive hard threshold on this 
value is meaningful only in the context of a theoretical 
model [l^ . we pragmatically choose a reference value of 
(j)T = 0.6 to distinguish a network with a strong small- 
world propensity from a network with weak small-world 
propensity iFig. 

The SWP has the critical ability to discern differ¬ 
ences in small-world structure across network densities 
(Fig. [IP). As density increases, the range of values 
spanned by the clustering coefficient and path length de¬ 
creases (Supplementary Fig. 1). The SWP is normalized 
by these ranges, minimizing the effects of network density 
on the calculation. The utility of this feature is particu¬ 
larly evident in comparison to the commonly used small- 
world index, tr, proposed by Humphries et al. [13)mi- As 
seen in Fig.[Tp-D, unlike the small-world index, the SWP 
retains a large dynamic range even as network density is 
increased, accurately pinpointing a small-world regime. 
(See the SI for a discussion of small-worldness in high 
density networks.) 

B. Generating Weighted Small-World Networks 


where 


and 



Ac 


Ciatt Cobs 
Ciatt Crand 


Ar = 


Lobs Lrand 
Liatt Lrand 


( 1 ) 

( 2 ) 

( 3 ) 


The ratio Ac /l represents the fractional deviation of the 
metric (Cobs or Lobs) from its respective null model (a 
lattice or random network). [13 See arrows in Fig. [TJ\ 
for a schematic of the deviation, and see Materials and 
Methods for mathematical definitions of C and L. 

In a standard Watts-Strogatz model, we observe that 
SWP is maximal for network configurations with the 
greatest small-world characteristics (Fig. [T]4-B). For 
small p, the SWP is low, driven by a high path length 
compared to that of a random network (resulting in a 
large Ac). However, as p increases, the path length 
quickly becomes closer to that of a random network, 
while the network retains a high clustering coefficient, 
similar to a lattice. Thus, when the SWP is maximal 
for p w 0.02, we see equal contributions from Ac and 
Ac. As p increases further, the network becomes in¬ 
creasingly random: the path length remains small, and 
the low SWP is now driven by the lack of local clustering 
(high Ac). 


To extend the SWP to weighted networks, it is nec¬ 
essary to define models of weighted lattice and random 
networks. Here, we take inspiration from the importance 
of physical distance in spatially embedded networks [H, 
and the observation that in many real physical and bi¬ 
ological networks, the strength of an edge is inversely 
correlated with the physical distance between nodes, i.e., 
nodes located near one another tend to be linked by 
stronger edges than nodes located far from one another 
[20l - l^ . By incorporating physical space, we define a re¬ 
wiring mechanism that allows a weighted network to be 
manipulated from a lattice network to a random network 
while maintaining the distribution of edge strengths: 

1. Begin with a network of N nodes that are arranged 
on a lattice and connected to all neighbors within 
a radius, r. 

2. Assign edge weights Wij according to the distance 
between nodes dij: 

niij = Umax dij (4) 

where Dmax = rnax{dij} -I- i5 is the maximum pos¬ 
sible distance between two nodes plus a single unit 
of measurement of the lattice spacing. [2^ ’ [^ 

3. Re-wire each edge with probability, p, retaining the 
weight of the edge. 

Following network formation, we assess small-world 
structure as a function of the rewiring probability by 
computing weighted versions of the clustering coefficient. 
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FIG. 1. Small-World Propensity in binary networks. 

(A) Clustering coefficient and path length as a function of 
the rewiring parameter, p, for a standard Watts-Strogatz for¬ 
mulation of a small-world network with N = 1000 nodes and 
r = 5. The dashed horizontal lines mark the baseline value 
of the clustering coefficient in a similar lattice network (top) 
and the baseline value for the path length in a similar random 
network (bottom). (B) SWP calculated for the same network 
as in panel (A). Error bars represent the standard error of the 
mean calculated over 50 simulations, and the shaded regions 
represent the range denoted as SW if using a threshold value 
of 4>t = 0.6. (C) SWP as a function of network density (in¬ 
creasing r for N = 1000 nodes). (D) Small-world index for 
the same networks as in panel (C). 

Cun and characteristic path length, (see Fig. 03 and 
Materials and Methods for definitions). These variables 
monitor a transition of our weighted networks through 
the small-world regime in a manner that is highly simi¬ 
lar to that observed in the original Watts-Strogatz pro¬ 
cedure for binary networks (compare Figs. [TJ\ & 03). 
We quantify this transition by computing a weighted 
version of the SWP using Ciatt = C^ip = 0) and 
Lrand = (p = 1). We observe that the weighted 
SWP behaves similarly to the unweighted SWP (com¬ 
pare Fig. [Tj3 & Fig. 0^), and maintains a large dynamic 
range as network density is increased (Fig. O^)- 


C. Mapping Real-World Observations to Theory 

To compute the SWP in a real-world network, it is 
necessary to determine the appropriate comparable lat¬ 
tice and random networks. To account for the effects of 
network density, we construct lattice and random net¬ 
works that maintain the network density - number of 


FIG. 2. Small-World Propensity in weighted networks. 

(A) Generation of weighted small-world networks. After 
building a lattice, the edges are weighted by distance such that 
close edges have a higher strength than distant edges. These 
edge weights are then retained as links are rewired with a 
probability, p, to create a weighted small-world network. (B) 
Weighted clustering coefficient and weighted path length as 
a function of the rewiring parameter, p, for a weighted for- 
mnlation of a small-world network with N = 1000 nodes and 
r = 5. (C) Weighted SWP calculated for the same network 
as in panel (B). Error bars represent the standard error of the 
mean calculated over 50 simulations, and the shaded regions 
represent the range denoted as small-world if nsing a thresh¬ 
old value of 4>t = 0.6. (D) Weighted SWP as a function of 
network density (increasing r for N = 1000 nodes). 


nodes and edge weight distribution - observed in the real- 
world network. Specifically, we construct a comparable 
weighted lattice by arranging the observed edge weights 
such that the edges that correspond to the smallest Eu¬ 
clidean distance between nodes are assigned the highest 
weights (Supplementary Fig. 2). For example, for a ID 
lattice with N nodes and unit spacing between nodes, we 
have N edges with a Euclidean distance of d = 1 between 
nodes. We therefore rank the observed edge weights by 
decreasing strength and randomly distribute the connec¬ 
tions with the N highest weights among the edges repre¬ 
senting d = 1. We proceed to distribute the next — 1 
edges of highest weight among the edges of the lattice cor¬ 
responding to d = 2, and we continue to proceed in this 
manner until the total number of edges in the real-world 
network have been placed in the lattice. To create a com¬ 
parable random network, the observed edge weights are 
randomly distributed among the N nodes of the network. 
These reference networks are then used to calculate the 
SWP of the observed network. 
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D. Synthetic Benchmark Networks 

To further assess the validity of the SWP, we examine 
its performance on weighted networks with well studied 
structure: fractal hierarchical (FH) and modular small- 
world (MSW) networks (see Materials and Methods). 
While both types of networks have a modular structure 
and therefore dense local clustering, the networks differ in 
terms of their path length: the MSW networks incorpo¬ 
rate random shortcuts of moderate connection strength 
between modules, decreasing the weighted path length, 
while the FH networks are composed of hierarchically or¬ 
ganized modules that are weakly interconnected, increas¬ 
ing the weighted path length (see Fig.[3]4). It follows that 
MSW networks have greater small-world structure than 
FH networks. 

In Fig. [21 we compare the ability of three metrics 
to detect small-world structure in these networks con¬ 
structed at various densities: (i) the weighted SWP com¬ 
puted on the true networks, and (ii) the un-weighted 
SWP and (iii) Humphries’ small-world index computed 
on binarized versions of the true networks. In the FH 
networks, we observe that both the weighted and un¬ 
weighted SWP values are relatively independent of den¬ 
sity, while the small-world index is strongly dependent 
on density. Moreover, the small-world index classifies low 
and medium density FH networks as having SW proper¬ 
ties (cr > 1), despite the long path lengths characteris¬ 
tic of these networks. The SWP instead reflects our ex¬ 
pectations, and remains below the small-world threshold 
{(j)T = 0.6) even in the case of binarized networks. Impor¬ 
tantly, when edge weights are retained, the SWP reveals 
that the FH networks are quite different than their MSW 
counterparts and do not show a small-world structure. 

In the MSW networks, we observe that the SWP values 
are sensitive to a meaningful change in the network struc¬ 
ture as density increases. Specifically, the path length 
becomes closer to that of a random network, but the ad¬ 
dition of shortcuts between modules reduces the cluster¬ 
ing coefficient thereby increasing the deviation from the 
clustering coefficient of a comparable lattice network Ac. 
This change is reflected by decreasing SWP. The small- 
world index is less sensitive to the deviation in clustering, 
and classifies all MSW networks as having small-world 
structure. 


E. Real-World Brain and Neuronal Networks 

The quintessential example of a small-world network 
in biology is the neuronal network of C. elegans In¬ 
deed, brain networks more broadly have been described 
as small-world for the past decade [1^ [2^. Yet, with 
advanced imaging and methodological techniques that 
provide more detailed measurements of these data with 
varying network densities and fine-scale estimates of edge 
strength, the assumptions of small-worldness have been 
brought into question Here we aim to resolve this 


A Low Density Medium Density High Density 



Low Med High Low Med High 


FIG. 3. Benchmark networks. (A) Adjacency matrices 
depicting network organization for increasingly dense Frac¬ 
tal Hierarchical (FH) and Modular Small-World (MSW) Net¬ 
works. (B) Comparison of the unweighted SWP, weighted 
SWP, and small-world index when applied to the matrices in 
panel (A). The dashed line represents the chosen threshold 
for indicating small-world structure {(j>T = 0.6 for SWP or 
or = 1 for the small-world index). Stars denote networks 
classified as small-world. 


controversy by applying the SWP to a representative set 
of structural and functional brain networks from several 
species whose properties have been well-studied in pre¬ 
vious literature, and whose edges are both binary and 
weighted (Fig. 01^). The three weighted networks are (i) 
a structural network representing the number of white 
matter tracts connecting 83 brain regions obtained from 
human diffusion spectrum imaging (DSI) data [l^, (ii) 
a functional network given by correlations between the 
blood oxygen level dependent signal of 638 brain regions 
measured using resting state functional magnetic reso¬ 
nance imaging (rs-fMRI) [2^, and (iii) a structural net¬ 
work of the cat cortex obtained from tract-tracing studies 
between 52 brain regions [s^. The two binary networks 
are (i) the neuronal network of C. elegans representing 
syna ptic and gap junction connections between 279 neu¬ 
rons [3l| , and (ii) a structural network derived from tract¬ 
tracing studies between 71 regions in the macaque cortex 


We observed that all networks displayed relatively 
large SWP values, indicating that brain and neuronal 
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FIG. 4. Real-world brain networks. (A) Adjacency ma¬ 
trices for three weighted brain networks (upper) and two bi¬ 
nary brain networks (lower). All matrices are symmetric. (B) 
SWP for the matrices shown in panel (A). The dashed line 
denotes the (pT = 0.6 threshold. (C) Breakdown of the SWP 
into the individual contributions from the clustering coeffi¬ 
cient and path length. High values of Ac or Ac indicate a 
large deviation from the comparable benchmark value, which 
results in a reduction of the SWP. 


networks in general do indeed display small-world prop¬ 
erties (Fig.0|3). Surprisingly, the network with the lowest 
SWP, sitting just below the (pr = 0.6 threshold, is the 
neuronal network of C. elegans, the quintessential exam¬ 
ple of a small-world network in biology @ . We explored 
this surprising result by examining the contributions of 
the clustering (Ac) and the path length (Ac) to the 
SWP value (Fig. HP). In contrast to the other brain net¬ 
works studied, the neuronal network of C. elegans dis¬ 
plays drastically different contributions: a remarkably 
high Ac indicating divergence from a lattice network, 
and an exceptionally low Ac indicating similarity to a 
random network. The observed SWP is therefore driven 
almost entirely by the short path length, indicating that, 
in reality, this canonical example of a small-world net¬ 
work, does not strongly embody small-world principles. 


III. DISCUSSION 

As petabytes of rich new data quickly file onto grow¬ 
ing storage farms, complex interconnected systems are 
being increasingly queried. Using these data to charac¬ 
terize and manipulate systems requires the development 
of network-based analysis tools that accurately account 
for the newly measured features of the data. Here we 
offer a novel network statistic (SWP) that accurately 
quantifies small-world structure in networked systems, 
that is agnostic to nuisance variables such as density 
and exquisitely sensitive to critical variables such as the 
strengths of connections between nodes. This enterprise 
has required the development of (i) a standardized pro¬ 
cedure for generating weighted small-world networks, (ii) 
a weighted extension of the SWP, and (iii) a stringent 
and generalizable method for mapping real-world data 
onto the theoretical model. We have illustrated the ap¬ 
plication and utility of these methods in the context of 
both benchmark and real-world data. The work repre¬ 
sents a single effort in the much larger space of meeting 
the demands of big data with increasingly sophisticated 
network-based tools. 


A. Small-World Statistics for Real-World Data 

The small-world index is the most common 

statistic to quantify small-world structure in binary net¬ 
works, but produces values greater than one for a large 
range of network topologies, obscuring accurate interpre¬ 
tations. Other measures monitor binary topologies with 
greater sensitivity but remain dependent on network den¬ 
sity and independent of edge weights [l^. More recent 
work has used weighted clustering coefficients and path 
lengths to compare observed networks to random net¬ 
works and to measure local small-world structure 
[S^ and topological dimension. However, the generaliz- 
ability of these methods has been hampered by the lack 
of a theoretical model to construct weighted small-world 
networks and to study the transition in and out of the 
small-world regime. (Although, see [s^ for a method 
to generate weighted scale-free small-world networks.) 
Our work directly fills this gap, allowing the creation 
of weighted small-world networks with varied degrees of 
SWP that can be used to compare and contrast net¬ 
work behavior and dynamics throughout the small-world 
regime. 

B. Binary Categories vs. Continuous Narration 

Recent network-focused efforts in the applied mathe¬ 
matics, physics, computer science, and engineering com¬ 
munities evidence the age-old tension to retain model 
simplicity while maximizing pragmatic utility. Network 
statistics are no exception: we often wish to obtain bi¬ 
nary categorizations rather than continuous descriptions. 
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For example, one might wish to emphatically state that 
a network does or does not display community structure, 
rich-club architecture, or small-worldness. Yet, arguably 
a more interesting and useful statement might assess gra¬ 


dations of these properties in real-world systems [3^ l3 


The true power of the SWP lies not in its ability to define 
small-world structure, but instead in its ability to quan¬ 
tify and compare the continuous degree of small-world 
structure between different networks. Our work there¬ 
fore complements ongoing efforts to extend traditionally 
categorical distinctions to continuous measurments in the 
context of core-periphery structure and community 
structure 1^. 


C. Pragmatic Utility in Neuronal and Brain 
Networks 


small-world propensity, predominantly due to a lack of 
local clustering. We speculate that several key biolog¬ 
ical feautres of this network may explain this surpris¬ 
ing contradiction to the historical literature. The wiring 
diagram of C. elegans is (i) a micro-scale network, (ii) 
represents neuron-to-neuron connections throughout the 
entire body of the organism rather than only the head, 
and (iii) is drawn from a comparatively symplistic or- 
gananism, evolutionarily speaking . These features of 
scale, physical extent, and evolutionary class may drive 
toplogical properties away from the small-world architec¬ 
ture to enable a different class of neural functions than 
those associated with the brains of higher-order animals. 
Future work must more stringently examine other canon¬ 
ical examples of small-world networks in diverse real- 
world systems, and thereby build a more accurate assess¬ 
ment of the role of topology in complex system function. 


Network statistics that are independent of density are 
critical for network comparisons in real-world settings. In 
the context of brain networks, the need for such statistics 
is underscored by the growing interest in examining neu¬ 
ronal networks across (i) time in development and normal 
aging [ 3 ^ 113, (ii) health and disease Idp . Id^. and (iii) 
different stages of neurodegeneration [^. In these and 
similar contexts, it becomes very difficult to determine 
if observed differences in small-world structure are sim¬ 
ply a function of the differences in network density, or 
if they represent a true form of topological reorganiza¬ 
tion. By minimizing the effects of network density in the 
computation of the SWP, we allow for a more direct com¬ 
parison of the topological network structure across time, 
and between groups. 

Network statistics that are highly sensitive to edge 
weights provide increased sensitivity to network func¬ 
tion. In the context of brain networks, the need for 
such statistics has only recently been actively appreci¬ 
ated, as a growing body of literature demonstrates that 
both healthy and diseased brain function is differentially 
driven by strong versus weak connections [3 1^ ■ Weak 
connections have traditionally been ignored because of 
commonly applied thresholding techniques , but have 
recently been identified as potential biomarkers in psy¬ 
chiatric pathologies [ij and as predictors of cognitive 
function and fluid intelligence In the future, it 

will be interesting to build on weighted network diagnos¬ 
tics like the SWP to provide novel quantifications of weak 
connectivity and its role in cognition. 


D. Surprising Biological Insights 


IV. MATERIALS AND METHODS 

All analysis was done using MATLAB (MathWorks) 
and code to compute the SWP in real-world networks 
can be accessed at http://www.seas.upenn.edu/^dsb/. 


A. Clustering coefficient 


To calculate the clustering coefficient, as in Q, we first 
calculate the local clustering coefficient for each node, Ci 
and then define the clustering coefficient, C, to be the 
average of the local coefficients 

c = ( 5 ) 


For binary networks, the local coefficient, Ci, of each 
node. Hi, is the fraction of closest neighbors of Ui that 
are also connected. Multiple extensions of Ci to weighted 
networks have been proposed: 

Onnela et al. [^ : 


'^*.0 — _ 2 ^^ ’^^{WijWjk'Wik) ^ , 

Barrat et al. 1^: 


( 6 ) 




Cz R = 


Wii + Wik 


ki{h “ 1) 7^ 2(wi) 




^ij^jk^iki ( 7 ) 


Zhang et al. 


The canonical example of a biological small-world net¬ 
work is the wiring diagram of C. elegans. However, the 
observation of small-world structure in this organism has 
been built on a simplification of the weighted wiring dia¬ 
gram to a binary graph. Using the weighted SWP, we 
observe that in fact, this network displays very little 


Y.WijWjkWik 
k k 

In these equations, Wij is the strength of a connection 
between nodes i and j, Wij = Wij/max{w), ki is the 
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number of edges connected to rii, and = 1 if a con¬ 
nection exists between nodes i and j (aij = 0 otherwise). 
For all results reported in the main manuscript, we used 
the definition given by Onnela et al. (Eqn. [S]), but sim¬ 
ilar results were obtained when using other definitions 
(see Supplementary Fig. 3). 


B. Path length 

The characteristic path length for a network is given 

by 

where for a binary network, dij is the shortest path be¬ 
tween nodes i and j. For weighted networks, as in we 
define the distance between two nodes to be dij = 1/ Wij . 


C. Synthetic benchmark networks 

The synthetic benchmark networks were created using 
the Brain Connectivity Toolbox with modifications 
to make networks symmetric and with weighting schemes 
as defined in [s^. All synthetic networks have N = 1024 
nodes. 

Fractal hierarchical networks were created from the 
makefractalCIJ function provided in the toolbox, with 
mxJv = 10, E = 2, and sz-cl = 5,6, or 7 for low, 
medium, and high density networks, respectively. This 
results in networks with 6, 5, or 4 hierarchical levels with 
a base module size of n = 32,64 or 128. At the lowest 
level of the hierarchy, modules are fully connected, and 
connections are placed within each hierarchical level with 
a probability p = 2~\ where I is the hierarchical level. 


The weight of each connection Wij = pij, such that the 
weight of a connection is equivalent to the probability 
that a connection exists. Because this method creates a 
directed network, the matrices were symmetrized by se¬ 
lecting the upper triangle of the resultant matrix and us¬ 
ing these connections to create an undirected, symmetric 
network. This procedure resulted in final networks with 
weighted (binary) densities of 4.5(10.8)%, 9.1(18.7)%, 
and 17.9(31.3)% for low, medium, and high density FH 
networks. 

Modular small-world networks were created from the 
makeevenClJ function provided in the toolbox, with 
N = 1024, K = 65000,100000, or 150000, and sz.cl = 6 
for low, medium, and high density networks, respectively. 
This creates networks with 16 fully connected modules of 
size n = 64 and E = K — 64512 randomly distributed 
edges between modules. Connection strengths within 
modules were set to Wij = 1 while strengths of inter¬ 
module edges were set to Wij = 0.5. As with the FH 
networks, this method creates a directed network, so the 
matrices were symmetrized by selecting the upper tri¬ 
angle of the resultant matrix and using these connec¬ 
tions to create an undirected, symmetric network. This 
procedure resulted in final networks with weighted (bi¬ 
nary) densities of 6.1(6.2)%, 7.9(9.5)%, and 10.3(14.3)% 
for low, medium, and high density MSW networks. 
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I. EFFECTS OF DENSITY ON CLUSTERING AND PATH LENGTH 


As discussed in the main manuscript and shown in Supplementary Figure 1 for a Watts- 
Strogatz network, the range of values spanned by the path length and clustering coefficient 
decreases with increasing network density. Here we dehne the range of clustering or path 
length to be the difference between its value in lattice and random networks, 1 — C{p = 
1)/U(p = 0) or 1 — L[p = 1)/L{p = 0). In the main manuscript, we present results showing 
calculations of the SWP and small-world index for network densities in the range 0 — 20%. 
At a network density of 20%, we can see that the path length only spans approximately 30% 
of it’s original range. Furthermore, by the point at which the density increases to 40%, the 
path length is the same in lattice and random networks. This implies that the concept of a 
Watts-Strogatz small-world network is ill-dehned for these networks of higher densities, and 
caution should be used when trying to impose this formalism on high density networks. The 
effect of network density on path length and clustering further increases the need to take 
edge weights into account when quantifying network structure. If a network with many weak 
connections is binarized, the density of the network can become quite high which influences 
measurements of network properties. 


H. MAPPING REAL-WORLD DATA TO THE THEORETICAL MODEL 

In order to quantify small-world structure in real-world networks, it is necessary to de¬ 
hne a method for mapping the observed data to the theoretical model used to generate 
small world networks. Specihcally, the calculation of the SWP relies on the generation of 
a comparable lattice network (p = 0 in the theoretical models) and a comparable random 
network (p = 1 in the theoretical models). To control for network density when generating 
these comparable null models, we preserve the number of nodes and the distribution of of 
edge weights. As described in the main text, to construct a comparable weighted lattice, 
we build a ID network such that the edges that correspond to the smallest Euclidean dis¬ 
tance between nodes are assigned the highest weights, whereas to construct a comparable 
random network, the observed edge weights are randomly distributed among the nodes. In 
Supplementary Fig. |2l we give an example of the resulting comparable lattice and random 
networks for the Human DSI data set analyzed in the main manuscript. In the comparable 
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lattice network, the connections with the highest strength have been distributed along the 
diagonal of the adjacency matrix, while in the comparable random network, connections are 
randomly assigned throughout the matrix. In both cases, we ensure that the resulting ad¬ 
jacency matrix is free of self-connections and remains symmetric (reflecting the undirected 
nature of our initial network). 


III. SWP USING ALTERNATIVE METHODS OF WEIGHTED CLUSTERING 

Multiple methods of computing a weighted clustering coefficient have been proposed [ij- 
Q], and in the main manuscript we present results using the weighted clustering coefficient 
as dehned by Onnela et al. [l| (see Matericals and Methods). This particular algorithm for 
computing the weighted clustering coefficient reflects subgraph intensity and has the advan¬ 
tage of being computationally efficient. However, other dehnitions of a weighted clustering 
coefficient exist and reflect other features of weighted clustering that might be desirable in 
certain data sets jd]. Here, we show the transition in and out of the small-world regime for 
the same weighted small-world network depicted in Fig. 2 of the main manuscript, along with 
the corresponding SWP, calculated for the weighted clustering coefficient dehned by Barrat 
et al. (Supplementary Fig. |3j\-B) and Zhang et al. (Supplementary Fig. [3p-D). The 
resultant weighted small-world network displays a similar transition through the small-world 
regime when measured using these alternative weighted clustering measures, and the asso¬ 
ciated SWP is remarkably similar to that of the SWP obtained using the Onnela measure 
(Fig. 2B-C). 


[1] J.-P. Onnela, J. Saramaki, J. Kertesz, and K. Kaski, Physical Review E 71 , 065103 (2005). 

[2] A. Barrat, M. Barthelemy, R. Pastor-Satorras, and A. Vespignani, Proceedings of the National 
Academy of Sciences 101 , 3747 (2004). 

[3] B. Zhang and S. Horvath, Statistical Applications in Genetics and Molecular Biology 4 (2005). 

[4] J. Saramaki, M. Kivela, J.-P. Onnela, K. Kaski, and J. Kertesz, Physical Review E 75 , 027105 
(2007). 


3 





1.2 



FIG. 1. Effects of density on clustering and path length The range of the clustering co¬ 
efficient (blue) and path length (red) corresponding to a Watts-Strogatz network for increasing 
densities (increasing r). Error bars are shown in gray and represent the standard error of the mean 
calculated over 50 simulations. 



FIG. 2. Mapping real-world data to the theoretical model Adjacency matrices representing 
the mapping of the Human DSI data from the main manuscript (middle) to a comparable lattice 
(left) and random (right) network. In both the comparable networks, the original edge distribution 
is maintained, but connections are redistributed to create a weighted lattice or random network. 
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FIG. 3. SWP using alternative methods of weighted clustering (A) Weighted clustering 
coefficient as defined by Barrat et al. and weighted path length as a function of the rewiring 
parameter, p, for a weighted formulation of a SWN with N = 1000 nodes and r = 5. (B) Weighted 
SWP calculated using the clustering in (A) . (C) Weighted clustering coefficient as defined by 
Zhang et al. and weighted path length as a function of the rewiring parameter, p, for a weighted 
formulation of a SWN with N = 1000 nodes and r = 5. (D) Weighted SWP calculated using 
the clustering in (C) . Error bars represent the standard error of the mean calculated over 50 
simulations, and the shaded regions represent the range denoted as SW if using a threshold value 
of (l)T = 0.6. (A) and (C) are comparable to Fig. 2B in the main manuscript and (C) and (D) are 
comparable to Fig. 2C. 
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